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We study the conversion of the latent heat of a supercooled quark-gluon 
plasma into excitations of the pion field. Supercooling is predicted to occur 

P 



in the late stages of the evolution of a quark-gluon plasma produced in ener- 



getic heavy-ion collisions, if the QCD phase transition is of first order. The 

H : 

supercooled, chirally symmetric state, which contains potential energy asso- 
ciated with an energetically unfavorable field configuration, eventually rolls 
down to the true minimum of the effective chiral potential. When this motion 
is described in terms of the linear sigma-model, we find that the energy of 
the coherent a— field is very efficiently converted into pionic excitations due 
to anharmonic oscillations around the minimum. The system is expected to 

partially thermalize before its disintegration. 
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The process of hadronization being of essentially nonperturbative character is not well 
understood at present. This leads to large theoretical uncertainties in describing the late 
stages of temporal evolution of the quark-gluon plasma predicted to be produced in very en- 
ergetic heavy ion collisions. Recently, the hadronization of an extended quark-gluon plasma 
has been studied within the nucleation model [|I]-§J by considering the growth of hadronic 
bubbles in the quark-gluon plasma undergoing a first-order phase transition. Thermal equi- 
librium between both phases was assumed. The authors of refs. |l|J^] found that the conver- 
sion of the quark-gluon plasma into hadrons is slow with a characteristic time much greater 
than 10 fm/c, resulting in significant supercooling of the plasma phase. According to these 
model calculations the system may cool to T ~ 0.7T C before hadronizing, corresponding to 
a reduction by a factor three in the entropy density. Of course, the total entropy remains 
(at least) constant; the falling entropy density is compensated by spatial expansion while 
the system cools below T c . 

When the transition to the hadronic phase finally occurs it is, in general, accompanied by 
significant reheating due to the release of the latent heat associated with the first-order phase 
transition. Additional entropy may be created in this process, but not necessarily so [[|. If 
the chirally symmetric phase is supercooled to just such an extent that its entropy content 
at the temperature Tq < T c exactly equals the entropy content of the equilibrated hadronic 
phase at the higher temperature Tjj reached after reheating, no entropy is generated. If the 
symmetric phase is supercooled further, then additional entropy is produced. 

The present paper is an attempt to analyze the dynamics of this reheating process within 
the context of an effective Lagrangian for the chiral phase transition. It is inspired by the 
recent work by Kofman, Linde, and Starobinsky [Q , who studied the reheating after inflation 
in the early Universe. We here discuss a nonequilibrium scenario, where the hadronization 
coincides with the chiral phase transition. Our picture is closely related to that of Rajagopal 
and Wilczek who studied the excitation of long-wavelength modes during the decay of 
a supercooled, false chiral vacuum state. However, we are here concerned not with the 
formation of coherent field domains, which has been numerically studied extensively ^-0, 
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but rather with the excitation of quasi-thermal modes of the pion field. 

When the quark-gluon plasma is supercooled, the system is in a moderately excited state 
over a chirally symmetric false vacuum. It does not matter for our purposes whether the 
entropy carrying modes in this state are best described by nearly massless quarks or by 
hadrons which are light due to chiral symmetry restoration, as some have suggested [Q. 
What is important is that the number of light modes is much larger in the supercooled, 
chirally symmetric state than in the hadronic state with broken chiral symmetry into which 
it decays. 

Once the instability is triggered, the system rolls down to the true chirally asymmetric 
minimum of the effective potential. Since the oscillations around this minimum are an- 
harmonic, the coherent field is converted into pionic excitations. As we argue below, the 
conversion mechanism is fast and strongly couples the long wavelength modes to those of 
short wavelength. The pionic excitations created by this mechanism are expected to (par- 
tially) thermalize before freeze-out. However, we emphasize that these pionic excitations 
do not necessarily increase the total entropy of the system significantly. In fact, our model 
makes no statement at all about an entropy increase. Because we do not explicitly con- 
sider the process of decoherence of the excited pionic modes they carry, strictly speaking, 
no entropy. It appears quite likely that during the decoherence process entropy is being 
transferred from other hadronic modes into the originally coherent excitations of the pion 
field. What our model does illuminate, however, is the mechanism by which the latent heat 
of the vacuum is converted into excitations of the pion field and on which time scale this 
happens. 

In the following, we analyze this hadronization mechanism within the linear sigma model 
defined by the effective Lagrangian 

C{t, X ) = \d^ a d^a ~ \Ma ~ V 2 ) 2 + mlvX (1) 

where <p a = (x, 7?) and the parameters are A ~ 20, v ~ 90 MeV, m v m 140 MeV. It is 
important to stress that in the deconfined, chirally symmetric phase the fields cj) a are treated 
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as an effective description of quarks and gluons. The c/> a — field is thought to dominate the 
low-energy dynamics of quarks and gluons near the transition, see e.g. ||. Below the critical 
temperature the chiral symmetry is broken and one deals with the physical fields a = v o — x 
and tt, which correspond to the sigma mesons and pions, respectively. Here vo ~ v is the 
minimum of the potential in ([I]). 

When the system is significantly supercooled in the symmetric chiral phase, where a ~ v 
and tt pa 0, the energy density is dominated by the potential energy 

A 4 

£pot - ^ ■ 

If this energy is converted into thermal energy of the pion gas, i.e. 

Spot — tth — Yq j 

where the pion mass is neg lected,[] one finds the temperature T ~ 135 MeV, which is close 
to the freeze-out temperature of pions observed in relativistic nucleus-nucleus collisions 0. 
Thus, even the complete supercooling of the system with subsequent reheating is not in 
conflict with the data. Let us then discuss the scenario in more detail. 
The sigma and pion fields satisfy the following equations of motion: 

d 2 + 2\v 2 - 3\va + \a 2 + At? 2 ] a = Awvf 2 , (2) 

d 2 + m 2 - 2Xv a + Act 2 + A^f 2 ] tt = . (3) 

We use these equations to discuss the temporal evolution of the system which initially is in 
the chirally symmetric phase. The fields are treated as classical. Since the system is assumed 
to be significantly supercooled, the initial values of the fields are taken as a ~ v and tt pa 0. 
Then, the system appears on the top of a "Mexican hat" potential. When the a— field is 



1 Since the average thermal energy of a massless pion is 3T, the zero mass or ultrarelativistic 
approximation is not bad even for temperatures of the order of the pion mass. 
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rolling down to the potential minimum at a = 0, its amplitude decreases. Therefore one 
expects that v 2 a > va 2 > cr 3 , and as a zeroth approximation, we neglect in eqs. (|2], 0) the 
terms which are quadratic and cubic in the fields. In this way we get the equations 



d 2 + ml]a^ = 0, Id 2 + ml] 7r(°> = , 



(4) 



with m a = y2Xv = 600 MeV, which have plane wave solutions. Keeping in mind the initial 
conditions, we choose the solutions of eqs. (§]) as 



= a"ocos(m (T t + <p) 



7f <°> = 



(5) 



where the field a® is assumed to be homogeneous. The solution fl5|) describes undamped 
collective oscillations of the a— field around its minimum. Now we substitute the solutions 
(H) into eqs. (0, |^) and keep the terms no more than quadratic in the cr— field and linear in 
the 7r— field. Changing the time variable (2z = mat + tp) and writing down the equations 
for modes labeled by the momentum k we get 



d 2 

— + A a - 2q a cos(2z) 







(6) 



where 



and 



d 2 

dz 



— + A w - 2g 7r cos(2z) 



ni 1] (z) 



mi 



A = 4 m - + k2 

mi 



q n = 2— . 

v 



(7) 



(8) 



(9) 



Eqs. (§,0) are written down in the form of the well known Mathieu equation JT0[, which 
corresponds to the wave equation in elliptic coordinates. As discussed above o"o < v, thus 
< q a < 3 and < q w < 2, while A a > 4 and A w > 0.2 . 
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The solutions of Mathieu's equation are of the form 

F v {z) = e ivz P{z) , 

where the constant v, which is called the characteristic exponent, depends on A and q and 
P(z) is a periodic function with the period n. It is easy to see that for an infinitesimally small 
q there are narrow regions of resonant solutions of the equations (|6]j7]) for A pa N 2 where 
N is an arbitrary integer number. Then, the characteristic frequency of the wave equation 
is a multiple of the frequency of the stimulating force. Let us consider for illustration the 
first resonant region around A = 1. The parameters A, q and v are then connected by the 
equation ]nj 



For A = 1 one immediately finds 
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Decomposing v into the real and imaginary parts with |Re v\ 3> |Im u\, we get 

i 

Ira u pa ± — = q . 

As expected, we have found an unstable mode which generally grows exponentially in time. 

When the amplitude of the stimulating force represented by q increases, the narrow 
resonances change into wide resonant bands in the A — q plane. The imaginary part of 
the characteristic exponent v is negative there and the solutions of Mathieu's equation are 
unstable. For example, when q = 1 the resonant regions extend for A from the intervals (0, 
1.9), (3.9, 4.4), (9.06, 9.08) etc. As seen, the first two resonances are rather broad while the 
third and higher ones remain narrow for q = 1. 

Although our zeroth order solution (|^) is homogeneous, i.e. k = 0, it effectively couples 
(via eqs. (|6], [?])) not only to the long wavelength pion modes but also to those with |k| pa m a . 
Let us consider how fast these modes grow. For a given value of A and q, the characteristic 
exponents can be read from the charts presented in |ItJ. For q a = 3 and q n = 2, which 



are the maximal values, one finds for A a pa A n pa 4 the imaginary characteristic exponents 



equal 0.35 and 0.17, respectively. Then, the particle numbers proportional to | cr^ 1 ^ | 2 and 
|^W|2 g row j n ti me as g*/" 1- with r a = 0.9 fm/c and = 1.9 fm/c, where r _1 = m a lm v. 
The characteristic time r n of growth of the occupation of single-particle pion modes appears 
much smaller than the nucleation time found in [|I|. Therefore, the mechanism described 
here can indeed be responsible for hadronization of the quark-gluon plasma. 

The hadronization which proceeds simultaneously with the transition from the chirally 
symmetric to asymmetric state has been discussed in |T2j . The hadronization time has been 
identified there with the characteristic time of rolling down from the top of the "Mexican hat" 
to the potential minimum, which is of order 0.5 fm/c. In our opinion, we deal with "physical" 
pions only when the coherence of the a and tt fields breaks down due to the anharmonic 
oscillations around the true vacuum state. Therefore, we identify the hadronization time 
with TV found above which equals about 2 fm/c. 

An important feature of our hadronization scenario is the strong coupling of the soft 
coherent modes to those with |k| « m a . This implies that the spectrum of produced 
pions is rather broad, similar to a thermal spectrum at T « m n . Therefore, we expect at 
least partial thermalization of the system. The reheating or even overheating due to the 
rapid release of the latent heat has also been advocated [13| within the nucleation model 
|Ij] . The thermal spectrum of hadrons observed in relativistic heavy- ion collisions supports 
such an expectation. It might be of critical importance for the so-called Disoriented Chiral 
Condensates actively discussed recently f|§,|nj • 

It has been suggested that coherent domains of the pion field can appear in the nonequi- 
librium chiral phase transition of the quark-gluon system created in high energy collisions. 
The phenomenon, analogous to the formation of misaligned domains in a ferro magnet, would 
be observed by the coherent pion emission when the domains relax to the ground state with 
vanishing pion field. Specifically, one expects significant fluctuations in the number of neu- 
tral and charged pions. The soft pions are expected to be particularly sensitive to the domain 
formation. Numerical simulations of the linear sigma model represented by the Lagrangian 
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(H) have shown |5]|| that the disoriented condensates can indeed appear under favorable 
conditions. 

As the reheating proceeds in the hadronizing system, it will generate additional back- 
ground on which it is more difficult to observe the pions from coherent domains. We note 
that some of the numerical simulations do not fully take into account the short wave- 
length modes excited by eqs. (f|[7|) due to the finite lattice spacing, typically of 1 fm. In the 
work of Asakawa et al. @ where the lattice spacing is 0.25 fm, there are indications that 
reheating occurs. 

In summary, if the quark-gluon plasma produced in relativistic heavy-ion collisions is 
significantly supercooled at the late stage of its evolution, its thermal energy is converted 
into potential energy of the a— field. The system then rolls down from the chirally symmetric 
to the asymmetric state and the hadronization coincides with the chiral phase transition. The 
physical pions emerge when the system anharmonically oscillates around the true vacuum. 
Due to the efficient resonant coupling between the soft and hard modes, the system is 
expected to thermalize again, at least partially, on a short time scale. 
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